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BLOW-UP FOR REALIZING HOMOTOPY CLASSES IN THE 
THREE-BODY PROBLEM 


Abstract. This expository note describes McGehee blow-up m in its role 
as one of the main tools in my recent proof with Rick Moeckel m that every 
free homotopy class for the planar three-body problem can be realized by a 
periodic solution. The main novelty is my use of energy-balance to motivate 
the transformation of McGehee. Another novelty is an explicit description of 
the blown-up reduced phase space for the planar N-body problem, N > 3 as 
a complex vector bundle over [0, oo) X CP J ' V where r G [0, oo) measures the 
size of a labelled planar N-gon and [s] E CP 2 ' 2-2 describes its shape. 


1. Introduction 

Deleting collisions endows the configuration space of the planar three-body prob¬ 
lem with non-trivial topology. Modulo rotations, this space is homotopic to a two- 
sphere minus three points and so has a large set of free homotopy classes of loops 
We call these classes the “relative free homotopy classes. We say that a solution is 
“relatively periodic" if it is periodic modulo rotation, or equivalently, if it is periodic 
in some rotating frame. 

Theorem 1 ((i?M) 2 2~ll). For equal or near-equal masses, and angular momenta 
J sufficiently small but nonzero, every relative free homotopy class for the planar 
three-body problem is realized by a relatively periodic orbit for the Newtonian planar 
three-body problem having energy —1 and angular momentum J. 

Remark. “Relative free homotopy classes” are encoded by “reduced syzygy 
sequences” : periodic lists of l’s 2’s and 3’s such as 123232... where the symbol i 
indicates that the three bodies have become instantaneously collinear with mass i 
in between j and k with i,j,k a permutation of the symbols 1,2,3. See for example 
(f2T1) for details. 

History of the THEOREM. A well-known theorem in Riemannian geometry 
asserts that on a compact Riemannian manifold every free homotopy class of loops is 
realized by a periodic geodesic. Inspired by this basic geometric fact Wu-yi-Hsiang, 
in 1996, posed the question: “Is every free homotopy class realized by a (relatively) 
periodic solution to the planar Newtonian three-body equation?” For 17 years I 
tried to use variational methods to prove that the answer is “yes” when the total 
angular momentum J is zero. Finally, at the urging of Carles Simo, in October of 
2014,1 gave up on variational methods and tested the waters of dynamical methods. 
Almost as soon as I gave up I realized that Rick Moeckel had come within epsilon 
of proving theorem 1 in the 1980s [(CZI. El. and E3)]- 

My purpose in this note is to give an exposition of one of the principal tools in 
our proof, McGehee blow-up [T3], and a sense of how we use this tool to prove the 
theorem. The main novelty is the use of energy-balance to motivate the mysterious 
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transformation of McGehee. Another novelty is an explicit description of the blown- 
up reduced phase space for the planar N-body problem, N > 3. For further reading 
on McGehee blow-up we recommend Moeckel ES9, Pages 222-225 and Chenciner 0. 


2. Background: Equations and Solutions. 


2.1. The Equations. The classical three-body problem demands that we solve 
the system of non-linear ODEs: 

TOigi = F 2 1 + -F 3 i 

(1) W2<72 = Fi2 + F 3 2 

m 3 q 3 = -f 7 2 3 +F 1 i3- 

where 

q a ~ qb 


( 2 ) 


F^ = Gm a m.b- 


' ab 


is the force exerted by mass m a on mass mb and 

Tab = | q a -qb\,q a € and m a ,G > 0. 


Here a, 6 = 1, 2,3 label the bodies. The dimension d for us will be 2. (The standard 
value is d = 3. ) The m a represent the values of point masses whose instantaneous 

i2 

positions are q a (t). The double dots indicate two time derivatives: q = The 
constant G is Newton’s gravitational constant and is physically needed to make 
dimensions match up. Being mathematicians, we can and do set G = 1. 


2.2. The Solutions of Euler and Lagrange. The only solutions to the three- 
body problem for which we have explicit formulae were found by Euler jlOj and 
Lagrange [T2] in the last half of the 18th century. See figures [lj [2] Their solutions 
are central to our story. 

For Lagrange’s solution, place the three masses at the vertices of an equilateral 
triangle and drop them: let them go from rest. They shrink homothetically shrink 
towards their common center of mass, remaining equilateral at each instant. The 
solution ends in finite time in triple collision. This motion forms half of Lagrange’s 
triple collision solution. To obtain the other half of Lagrange’s solution use time- 
reversal invariance to continue this solution backwards in time. In the full solution 
the three masses explode out of triple collision, reach a maximum size at the instant 
at which we dropped the three masses, and then shrink back to triple collision, 
staying equilateral throughout. A surprise is that the Lagrange solution works 
regardless of the mass ratios m 3 : m 2 : m 3 . 

For Euler’s solutions, place the masses on the line in a certain order: qi < qj < q^ 
so as to form a special ratio qu — qj : qj — qi . (This special ratio depends on the mass 
ratios and also the choice of mass mj on the middle and is the root of a fifth degree 
polynomial whose coefficients depend on the masses.) Again drop them. They stay 
on the line as they evolve and again the similarity class of the (degenerate) triangle 
stays constant: this ratio of side lengths stays constant. ( In case the two masses 
at the ends are equal then the special ratio is 1 : 1: place mj at the midpoint of mi 
and TOfc. ) 

The solutions just described are part of a family of explicit solutions discovered 
by Euler and Lagrange. For every one of the solutions in these families the similarity 
class formed by the three masses stays constant in time during the evolution. Each 
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mass moves on its own Keplerian conic with the center of mass of the triple as 
focus, the solutions described i above being the special case of degenerate (colinear) 
ellipses. We derive these families analytically in section [4. 3. 1| below. 



Figure 1. A Lagrange Solution. 


All together these solutions form five familes. The corresponding shapes are 
called “central configurations”. The Lagrange solutions count as two, one shape 
for each orientation of a labelled equilateral triangle. The Euler solutions count as 
three, one for each choice of mass in the middle. 

For almost all (Newtonian) time the solutions of theorem 1 are very close to one 
of the three Euler solutions. The Lagrange solutions act as bridges between various 
Eulers. 
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Figure 2. An Euler Solution. 

3. Shape sphere. Blow-up and reduction, first pass. 

A basic aid to understanding the planar three-body problem is the shape sphere , 
a two-sphere whose points represent oriented similarity classes of triangles. At 
each instant of time three bodies form the vertices of a triangle. Call two triangles 
“oriented similar” if one can be brought to the other by a composition of transla¬ 
tions, rotations, and scalings. The resulting space of equivalence classes forms the 
shape sphere. See figure [3] This sphere has 8 marked points, the 5 central con¬ 
figurations just described L + ,L_,Ei,E 2 ,E 3 and 3 binary collision points labelled 
B 12 , B 23 , B 3 i . The sphere’s equator represents the space of collinear triangles. 
The 3 binary collision points, and 3 Euler central configurations lie on this equator, 
interleaved so as to be alternating. 

The earliest occuring picture of the shape sphere in the context of celestial me¬ 
chanics with which I am familiar is m . You will find a detailed exposition of the 
shape sphere and its relation to the three-body problem in [22l . 

We now summarize how the shape sphere arises out of the three-body problem. 
The configuration space for the 3 body problem, with collisions allowed, is C 3 with 
q = ( 91 , 92 , 93 ) £ C 3 representing the 3 vertices of the triangle - the positions of 
the 3 bodies. We have identified C with R 2 in the standard way: x + iy £ C 
corresponds (x,y) £ R 2 . A standard trick from Freshman physics allows us to 
restrict the problem to the center-of-mass zero subspace: 

E cm = {q £ C 2 : TO 191 + m 2 q 2 + m 3 q 3 = 0} = C 2 C C 3 . 

(See the beginning of section [ 5 ] below.) In E cm the binary collision locus become 
three complex lines which intersect at the origin 0. The origin represents triple 
collision. The masses endow C 3 with a canonical metric called the “mass-metric” 
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Figure 3. The shape sphere. Lagrange points, Euler points, and 
collision points marked. The equator consists of collinear triangles, 
figure courtesy of Rick Moeckel 


(eq. 0 ) and relative to that metric the distance from triple colliison is given by r 
where 

r 2 = mi|gi | 2 +m 2 |g 2 | 2 + m 3 |g 3 | 2 . 

(See eq. Q.) Take the sphere 

{r = 1} : = S 3 C C 2 = E cm . 

Because the three-body equations are invariant under rotations they descend to 
ODEs on the quotient of C 2 = E cm by the group S 1 of rotations. This quotient 
space is topologically an R 3 . To understand this quotient note that the rotation 
action leaves r unchanged but moves points on S 3 around according to {Z\ 1 Z 2 ) i—»• 
[uZ i, uZ 2 ), u £ S 1 C C. (Here Z i, Z 2 are any complex linear coordinates for E cm .) 
This is the circle action used to form the Hopf fibration: 

Hopf : S 3 —» S 3 /S 1 = S 2 = shape sphere . 

Points of the quotient R 3 represent oriented congruence classes of triangles: planar 
triangles modulo translation and rotation, but not scaling. Express R 3 in spherical 
coordinates (r,s),s € S 2 . Then the origin r = 0 corresponds to triple collision. 
A point s on the sphere represents a ray rs, r > 0 of triangles all having the 
same shape. The collision locus C = {ri 2 = 0 or r 23 = 0 or r 33 = 0} is then 
represented by the three rays corresponding to the three binary collision points 
B 12 , -B 23 , B 3 i G S 2 . 

Newton’s equations break down at triple collision r = 0. McGehee blow-up is a 
change of variables ( equations @) which converts Newton’s equations to a system 
of ODEs which is well-defined when r = 0. The locus r = 0 in the new variables 
is called “the collision manifold” and forms a bundle over the shape sphere. The 
blown-up system of ODEs has exactly 10 fixed points, all on the collision manifold 
, with a pair of fixed points lying over each of the five central configurations. For 



6 BLOW-UP FOR REALIZING HOMOTOPY CLASSES IN THE THREE-BODY PROBLEM 


a chosen central configuration, one element of the pair corresponds to the homo- 
tlietic arc incoming to triple collision, as in our original description of the Lagrange 
solution, while the other element of the pair corresponds to the initial segment of 
that solution which explodes out from triple collision. 

The 10 fixed points on the collision manifold have stable and unstable manifolds, 
parts of which stick out of the collision manifold, and which intersect in complicated 
ways, as per the Smale Horseshoe and heteroclinic tangles. See figure [5] Moeckel 
investigated these manifolds and their relations in seminal works US], US], EU, 
m, nzi> and m where he proved existence of “topological heteroclinic tangles” 
between them. 

One finds the following abstract graph 



in several of these papers of Moeckel am p. 53, Theorem l 7 , and m whose 
Figure 2 becomes our graph after deleting the vertices labelled with Bs (for binary 
collision) as well edges incident to them). Moeckel’s theorem in [17], based on the 
intersections between stable and unstable manifolds of the 10 fixed points, asserts 
that all paths in this graph are “realized” by solutions to the three-body problem 
provided the angular momentum, energy and masses are as per theorem 1. Embed 
this graph in the shape sphere as indicated by figure [3] Call the embedded graph 
the “concrete connection graph”. 

The dynamical relevance of the concrete connection graph has to do with the 
Isosceles three-body problem. When two of the masses are equal, say mi = m 2 , 
then the isosceles triangles r 13 = ?~23 form an invariant submanifold of the three- 
body problem whose dynamics is called the “Isosceles three-body problem”. These 
Isosceles triangles form a great circle passing through both Lagrange points, the 
binary point B 23 , and the Euler point E\. If all three masses are equal we have 
have three Isosceles subproblems represented by three great circles on the shape 
sphere. Take one-half of each great circle, namely that half whose endpoints are 
the two Lagrange points and which contains the Euler point. In this way we form 
the concrete connection graph in which the edges are parts of the Isosceles great 
circles. 

Observe that the shape sphere minus the three binary collision points retracts 
onto the concrete connection graph. Theorem 1 follows immediately from this 
observation and Moeckel’s theorem refered to above, once we know that the realizing 
solutions of Moeckel’s theorem, projected onto the shape sphere, stay C'°-close to 
corresponding edges in the concrete connection graph. For a few more details see 
the final section of this article. 

4. Set-up. Blow-up 

It is no more work to perform the blow-up for the N body problem in d- 
dimensional Euclidean space, rather than our special case of the three-body problem 
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Figure 4. The concrete graph, embedded in the shape sphere, 
in the plane. The d-dimensional N-body equations are: 

(3) 'W'aQa = ^b^a-^ba : Qa £ 

with the forces Ft, a as above. 

4.1. Metric Reformulation. Let 

E = {R d ) N 

denote the N-body configuration space. Write points of E as q = (qi ,..., qn) and 
think of the points as the N-gons in d-space. The masses endow E with an inner 
product, called the mass inner product: 

(4) (q, v) = E m a q a ■ v a 
so that the standard kinetic energy is given by 

(5) K=^{q,q). 

Let V be the gradient associated to this metric: df q (v) = (V/(g),u), so that 
(V/) 0 = Then the N-body equations take the simple form 

(6) q = VU(q) 
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Figure 5. The equilibria arising upon blow-up and relations be¬ 
tween their stable and unstable manifolds. The purple and green 
arrows are the rest cycles described in figure pTTj 


where U is the negative of the standard potential V : 


(7) 


U = -V = T, 


a<b 


m a m b 

Tab 


the sum being over all distinct pairs a , b. As is well known, the total energy 
conserved (i.e constant) along solutions 


H = K-U = K + V. 


We use 

(8) r = y/{q,q) 

to measure the size of our configuration q = (qi ,..., gjy). Lagrange proved that 

r 2 = ’£ a <bm a m b rl b /T,m a 

provided we are in center of mass coordinates: T,m a q a = 0. Then 
r = 0 <t=> total collision: all masses coincide 


while 


U = oo <t=> some collision: some pair of masses coincide. 
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Exercise 1. Use the metric reformulation of Newton s equations eq the fact 
that U is homogeneous of degree —1 and Euler’s identity for homogeneous func¬ 
tions to derive the “virial identity”, also known as the Lagrange-Jacobi identity: 
d 2 {r 2 )/dt 2 = AH + 2U. Also show that 2 H + U = H + K. 

4.2. McGehee transformation via Energy Balance. The key property of the 
potential energy —U, as far as McGehee’s transformation is concerned, is that it is 
homogeneous of degree —1: U(Xq) = X~ 1 U(q), or 

q ha Xq =+> U i—^ X~ 1 U. 

Our guiding principle in deriving the McGehee blown-up equations is to require that 
the kinetic energy K scale the same as potential energy, so that the total energy 
balance under scaling and thus has a scaling law. We call this principle “energy 
balance”. Then K ha A _1 K and since K is quadratic in velocities v velocities must 
scale by 

v = q ha A -1 / 2 il 

How must time scale? Since dq ha Xdq and v = dq/dt, we see that for a power 
law scaling dt ha A a dt to yield v ha A~ 1//2 u we must have a = 3/2. Summarizing, 
our space-time scaling law must be 

(9) q ha Xq, dt ha A z t 2 dt, 
which induces the desired scalings 

v ha X~ 1/2 w,{U,K,H) ha A _1 (I7 ,K,H) 

Exercise 2. Show that q(t) solves if and only also q\(t) := Ag(A~ 3 / 2 t) solves 
|6|b Explain how the exponent —3/2 in this transformation-of-paths formula arises 
from the +3/2 in the time part of the scaling law of eq &■ 

McGehee’s genius was to rewrite Newton’s equations, as much as is possible, in 
scale invariant terms. We cannot completely get rid of scale, but we can encode 
scale in the single size variable r = (q, q) introduced earlier and through which 

we remove scale from the remaining variables: 

(10) q 

(11) V 

(12) dt 

These relations define the McGehee transformation ( q,w,t ) ha (r, s,y;r). Observe 

that s lies on the unit sphere r = 1 in the configuration space, 

sGS = S dAr_1 ={r = l}cE 

so that (r, s ) are spherical coordinates on E. We sometimes refer to s as the shape 
of the configuration q. 

Exercise 3. Write ' for -4- = r 3 / 2 ^-. Show that McGehee’s transformation trans¬ 
forms Newton’s equations Jhj) to the equations 

r' = rv 


= rs 



= r 3 ' 2 dr 


(13) 


s' = y — us 
li = VC/(s) + \vy 


where v = (s,y ). These equations are the McGehee blown-up equations. 



10 BLOW-UP FOR REALIZING HOMOTOPY CLASSES IN THE THREE-BODY PROBLEM 


In the last equation Vt/(s) £ E is the same gradient as in Newton’s equation Q, 
only that restricted to points s of the sphere {r = 1}. The blown-up equations are 
analytic and extend analytically to the total collision manifold r = 0. For N > 2 
the equations still have singularities due to partial collisions eg ri 2 = 0, at which 
Vf7(s) still blows up. 

Definition 1. The “extended collision manifold” is the locus r = 0 for the blown-up 
phase space [0, oo) x S x R. dN of McGehee. 

The first of the three blown-up ODEs asserts that the extended collision manifold 
is an invariant submanifold. On the extended collision manifold the flow is non¬ 
trivial , as a glance at the last two equations shows. Away from the extended 
collision manifold, the blown-up equations are equivalent to Newton’s equations. 
What have we gained by adding this collision manifold? 


4.3. Equilibria! The first thing one learns in a class in dynamical systems is to 
look for equilibria. But Newton’s equations have no equilibria! N stars cannot 
just sit there, still, in space. Adding the extended collision manifold thru blow-up 
introduces equilibria. When N = 3 these equilibria correspond to the solutions of 
Euler and Lagrange described above. In the general case the equilibria correspond 
to “central configurations”. (Proposition [l] below.) 


Finding the equilibria. From the first of the blow-up equations (13) we see that 
at an equilibrium must lie on the extended collision manifold r = 0 (consistent 
with wha we just said about “stars cannot just sit there”). Plugging the second 
equilibrium equation 0 = y — iss into the third equation of the blown-up equations 
(131 yields the “shape equation” 


VU(s) = -\v 2 s. 

Taking the inner product of both sides of the shape equation with s and using 
Euler’s identity for homogeneous functions yields 

U{s) = ^u 2 


or v = ± sj2U{s). Now the gradient of the function r 2 at a point s £ S is 2s so 
that we can rewrite the shape equation as 

Vt/(s) = cV(r 2 ) , c = -\v 2 = - l -U{s). 

Think of c as a Lagrange multiplier. We have proved that the shape s of an 
equilibrium configuration must be a “central configuration” where: 


Definition 2. A central configuration is a shape s £ S which is a critical point of 
U restricted to the sphere r = 1. 


Conversely, for each central configuration shape s cc £ S' we obtain an equilibrium 
point ( r,s,y) = (0 ,s cc ,y) by setting y = vs cc . with v = ±^2 U(s cc ).. We have 
established 


Proposition 1. Equilibria of the blown-up equation are in 2:1 correspondence with 
central configurations. This correspondence associates to a given central configura¬ 
tion s cc the two equilibria ( r,s,y ) = (0, s cc , vs cc ), with v = zt^/2U(s cc )s cc ). 
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There is another way to arrive at central configurations in keeping with our 
original discussion of the Euler and Lagrange solutions. Make the ansatz: 

(14) q(t) = X (t)s 

where A(t) is a time dependent scalar and s £ S is constant. 

Exercise 4. Show that the ansatz \lf\j satisfies Newton’s equations if and only if 
s is a central configuration and A (t) satisfies the “Kepler problem”: 

(15) A = -/xA/|A| 3 m =^i/ 2 = [/(s) 


All solutions of the one-dimensional Kepler problem end in collision: A = 0. This 
value of A in the ansatz corresponds to total collision. The ansatz (14) with real 
scalar A(t) yields the Lagrange and Euler solutions to the three body equations 
which we first described above by “dropping” bodies. 


4.3.1. The Euler and Lagrange family. Planar problems. Assume we are in the 
planar case so that d = 2. Identify M 2 with C: (x,y) — > x + iy so that E = C N 
and so that complex scalar multiplication of s = (si,..., sn) £ S C C N by A £ C 
corresponds to scaling the N-gon s by the factor |A| while rotating it by Arg(X). 

Exercise 5. Show that NU(Xq) = j^'S/U(q) 

Exercise 6. Use exercise [5] to show that exercise [7] also holds in the case of the 
planar N-body problem, with now X (t) £ C a complex scalar. 


The solutions of exercise [6] are motions in which all N bodies move “homograph- 
ically” : meaning by scaling and rotating. Since A (t) parameterizes a conic, each 
body moves along a homographic conic A (t)s a , a = 1,2 ...,1V. We now have, for 
each planar central configuration, a family of solutions parameterized by the com¬ 
plex solutions to eq. (15), and varying from total collision solutions when A(t) £ M 
to circular motions A (t) = e lult £ S 1 C C. For fixed energy h we can think of the 
parameter of the family as the angular momentum J discussed below, with J = 0 
being total collision and the maximum or minimum value of J being the circular 
motion. 


4.3.2. An open problem. The potential U is invariant under rotations and transla¬ 
tions. Consequently, the central configurations as we defined them are not isolated, 
but come in families. 

Is the set of central configurations, modulo rotatations and transla¬ 
tions, a finite set? This problem is attributed to Chazy . See Albouy-Cabral 
[2] and for perspective and a recent survey. 

What is known. Some History. N = 3: Euler and Lagrange had established 
the complete list of central configurations as described here. N = 4: Albouy [1] 
classified the central configurations in the case of 4 equal masses two centuries two 
decades and a few years after Euler and Lagrange . One of his main achievements 
was to show that in the equal mass case the 4-body central configurations all have a 
reflectional symmetry. Eleven years after Albouy’s work Hampton and Moeckel m 
proved that the central configurations are finite (less than 1856 XXX in number!) 
N = 5: In 2012 Albouy and Kaloshin [3] proved that for N = 5 and away from an 
algebraic surface in the parameter space RP+ of mass ratios, the number of central 
configurations is finite. In 1999 Roberts [23] constructed examples for N = 5 with 
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one of the five masses negative in which the set of central configurations is infinite, 
underlining the subtlety of the problem. 


4.4. Linear and angular momentum. Besides energy, the only known constants 
of the motion for the general N-body problem are the components of the linear 
momentum 

P = T,m a v a 

and the angular momentum 

(16) J(q,v) = T,m a q a Av a 

These momenta are intimately connected to the fact that the group G of rigid 
motions acts by symmetries of Newton’s equations. 

Exercise 7. is orthogonal to the G orbit thru q £ K if and only if P{v) = 0 

and J(q, v) = 0 


4.5. Center of mass frame. A well-know argument using Galilean symmetry and 
found in essentially any introductory physics text allows us to suppose that all our 
solutions satisfy P = 0 and 

Zm a q a = 0 

In this case we say that we are in “center of mass frame” and we set 
E cm = {q £ E : Zm a q a = 0} £* R d ( N ~ 1 ' > . 


The infinitesimal generators of the translation action are the vectors q = (c, c,..., c) = 
“cl",c £ R d . Now E cm is precisely the orthogonal complement to the subspace of 
vectors of the form cl. This space of vectors corresponds to the generators of the 
translation group, or alternatively to the space of all total collision configurations. 

It follows that E cm realizes the quotient of E by translations and that in E cm only 
one point represents total collision: the origin. 

We can go to center of mass frame before or after blow-up, the result is the 
same, namely the system of ODEs (13) (with poles on the partial collision locus) 
but restricted to the the subset of variables 


(17) 

where 


(r,s,y) 


[0, Oo) X S C m X E cm 


S cm = {q£ E cm : ( 3 , s) = 1} = S^-U-i. 


4.6. Energy-momentum level sets and the Standard Collision Manifold. 

Because energy and angular momentum are invariant as we flow according to New¬ 
ton, by fixing their values h and Jo we obtain invariant submanifolds of phase 
space: 

M int (h) = {H = h,r > 0} 

and 

M int (h,J 0 ) = {H = h,r > 0} 

Energy and angular momentum are not defined at r = 0 so we have excluded r = 0. 
Set 

(18) M{h) = Closur e(M int (h)), M(h, J 0 ) = Closur e(M int (h, J 0 )), 

the closure being within within the blown-up phase space. We will need to under¬ 
stand the boundaries of these spaces, which is their intersection with the extended 
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collision manifold r = 0; in other words we must understand how these invariant 
submanifolds approach the extended collision manifold {r = 0 }asr—>-0 . 

The following notation will be useful in this endeavor. 

Definition 3. [Notation] For F = F(q,v) a homogeneous function on E x E 
write F for the scale-invariant version of F achieved by multiplying F by r~ a 
where a is the degree of homeqeneity of F with respect to our weighted scaling. 
Thus: F{q,v) = r a F(s,y). 


According to “energy balance” both the potential energy, kinetic energy, and 
total energy are homogeneous of degree — 1. Thus 

U{s) = rU{q) 

where U is homogeneous of degree 0 and can be viewed as a function on the sphere 
hi cm • And 

(19) H = rH 


where 

H{s,y) = ^(y,y) - U(s ) = I<(y) -U(s) 

and K, U are homogeneous of degree 0. The angular momentum is homogeneous 
of degree 1/2 so that 

(20) J = r 1/2 J(s,y) 


where J is scale invariant and equals T,m a s a A y a . 

If follows immediately from eq (19) that 

d{M{h)) = {H = 0,r = 0}. 


while using in addition eq (20) we see that 


d{M(h , 0)) = {H = 0, J = 0, r = 0}. 

These are basic important submanifolds so we give them separate names. 
Definition 4. The full collision manifold is M 0 = {H = 0,r = 0}. 
Definition 5. The “standard collision manifold” is the locus 

C:={ r = H = j = 0}. 


Thus the extended collision manifold contains the full collision manifold Mq 
which in turn contains the standard collision manifold C. The equilibria all lie on 
C. Another reason for the importance of the standard collision manifold C is a 
theorem of Sundman. 


Theorem 2. (Sundman) If r —> 0 along an honest solution, then J = 0 for that 
solution and hence that solution tends to C as r —> 0. Moreover, the solution tends 
to the subset of equilibria within C. 

Here we are using the hopefully obvious 

Definition 6. An “honest solution” to the blown-up equations is a solution such 
that r > 0. 
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The honest solutions are just the reparameterizations of solutions to our original 
Newton’s equations according to the blown-up time. 

Remark. The standard collision manifold C is the space most authors refer to 
when they speak of the “collision manifold” for the N-body problem. Clienciner 
(see also [7]) argues that the standard collision manifold is the dilation quotient of 
the N-body phase space. 


Exercise 8. 


Use eq (13) to show that 
d 

dr 


H = vH. 


^Lj = -\ v j 

dr 2 

hold everywhere on the blown-up phase space. 


It follows from this exercise that {H = 0} and {J = 0} are invariant manifolds, 
as are M 0 and C. 


5. Quotient by Rotations. 


Newton’s equations and their McGehee blow-ups (eq 13) are invariant under 


the group G of rigid motions and so descend to ODEs on the quotient space of 
their phase spaces by G. Working on this quotient instead of the original helps our 
intuition enormously, especially in the case N = 3 and d = 2. We describe the 
quotient and some aspects of the quotient flow. 

The group G of rigid motions is the product of two subgroups, the translation 
group and the rotation group. We have already formed the quotient of phase 
space by translations when we went to center-of-mass frame, i.e. by restricting 
to s,y £ E cm . To form the remaining quotient by rotations it is much cleaner 
to restrict to the planar case d = 2. Henceforth we assume that we are working 
with the planar N-body problem, d = 2. We identify R 2 with C as before. Thus 
E = C N and E cm = C^' 1 . Represent rotations as unit complex scalars u £ S 1 C C 
acting on (q, v) £ E cm x E cm by (q, v ) (uq, uv) and on McGehee coordinates by 
(' r,s,y ) i-A (r,us,uy). 


Definition 7. The blown up reduced phase space in the planar case is the quotient of 
the blown-up center of mass phase space [0, oo) x S cm x E cm = [0, oo) x S 2N ^ 3 xC N ^ x 
by the group of rotations. Upon deleting the collision locus C we denote this quotient 
by 

V N = ([0, oo) x ( S 2N ~ 3 \ C) x C^" 1 )/^ 1 . 

Momentarily forget the velocities v or y, and the deletion of the collision locus C 
in trying to understand the quotient. The circle action sends a blown-up configura¬ 
tion (r, s) to (r, us), s £ S cm . So we need to understand the quotient of the sphere 
S cm = S 2N ~ 3 by this action of S 1 . It is well known that this quotient S cm /S l is 
isomorphic to the complex projective space CP W_2 := P(E cm ) with the projection 
map S cm —» Scm/S 1 being the Hopf hbration. Hence the quotient of the (r, s) by 
S' 1 yields [0,oo) x CP^ -2 . 

To better understand the meaning of points of CP^ -2 , work with q £ E cm 
instead of q = s £ S cm , insisting only that g / 0 and now allowing the scalar u to 
vary over the larger group C* D S 1 of all nonzero complex numbers. The resulting 
quotient is well-known to be (C w_1 \ {0})/C* = CP^ -2 . The action of u £ C* on 
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q £ \ {0} is precisely the action of rotating and scaling the (centered) N-gon 

(J- 

Definition 8. The projective space CP^ -2 just constructed is called shape space. 
Its points represent oriented similarity classes of planar N-gons. 

We have realized the configuration part of the quotient after blow-up as [0, oo) x 
CP^" 2 where CP^" 2 is the shape space. When N = 3 the shape space is the shape 
sphere described above. 

Collision locus. 

The condition that a configuration q = (q lt ... ,qj\r) represent a collision is that 
q a = qb for some a ^ b,l < a,b < N. This condition is complex linear when viewed 
in homogeneous coordinates [<?i, qi, ■ ■ ■, <y/v] and so defines a complex hyperplane, 
a CP^ -3 C CP W ~ 2 . There are are (^) pairs (a, b) and so we have to delete 
( J ^) hyperplanes from our shape space. The union of these hyperplanes, viewed 
projectively, is the collision locus: 

C = {[<?] = [gi> < 72 , • • •, 9iv] e CP W ~ 2 :q a = q b some a/t}. 

We use the same symbol for the collision locus before or after quotient. 

Accounting for velocities. In the last few paragraphs above we dropped the 
velocity y. The quotient map (r, s) i-> (r, [s]) from [0, oo) x S cm —> [0, oo) x CP W ^ 2 
expresses [0, oo) x S cm as a principal S 1 bundle over [0, oo) x CP^ -2 . 

Now include the velocity y. The quotient procedure with y included is precisely 
the procedure used to construct an associated vector bundle to a principal bundle. 
Realizing this, we see that the quotient Vn is a complex vector bundle over [0, oo) x 
(CP^ -2 \C) whose rank is N — 1 - the fiber being coordinatized by y £ E cm . What 
is this vector bundle? 

Proposition 2. 

V N = [0, oo) x T(CP n ~ 2 \C)xl 2 

as a vector bundle over [0, oo) x (CP W ~ 2 \ C). The final R 2 factor is coordinatized 
by (u, J) where v = (s,y) represents the time rate of change of size and where 
J = ( is,y ) is also equal to r~ l ! 2 J off of r = 0 where J is the usual total angu¬ 
lar momentum of the system. The fiber variable tangent to shape space CP W ~ 2 
represents “shape” velocity. 

In the case of N = 3 we have CP^ 2 = CP 1 = S 2 , the shape sphere previously 
discussed in section [5] Then 

Vs = [0, oo) x T(S 2 \C)xR 2 = [0,oo)xRx (S 2 \C)xR 3 

where 

C = {Bi2,B 2 3, B 3 i} 

is the set of three binary collision points. 

5.0.1. Velocity (Saari) decomposition. Passing thru a configuration q £ E cm we 
have two group-defined curves: the scalings A q 1 A £ R of q and the rotations uq , u £ 
S 1 of q. The tangent spaces to these curves are orthogonal, and together with the 
orthogonal complement of their span they define a geometric splitting of T q E cm = 
E cm 

(21) T q E cm = 

( 22 ) 


(scale) + (rotation) + (horizontal ) 

R q ® iRq ® {u : J(q, v) = 0, v{q, v) = 0} 
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where 


Definition 9. The horizontal space at q is the orthogonal complement (rel. the 
mass metric) of the sum of first two subspaces R q and iR.q, i.e it is the orthogonal 
complement to the C-span of q. 

Refer to exercise [7] and the definition of v to see why the horizontal space at q is, 
as described above, the zero locus of J(q,v ) and v(q,v). 

Unit vectors spanning the scale and rotation spaces are s and is. Consequently, 
if we take & v £ T q E cm and decompose it accordingly we get 

(23) v = (s, v)s + ( is , v)is + Vhor 
and the scale invariant version: 

(24) y = us + Jis + y hor ', u = (s,y), J = (is,y) 

where the subscript “hor” on v and y denote their orthogonal projections onto the 
horizontal subspace. 

Remark D. Saari pointed out the importance to celestial mechanics of the 
horizontal-vertical splitting of eq ( [23] ) and hence this splitting is often called the 
“Saari decomposition”. 


5.0.2. Proof of proposition^ The decomposition (eq (24)) of y is 5' 1 -equivariant. 
The coefficients of the first two terms v and J = ( is 7 y ) are S^-invariant functions 
and so are well defined functions on the quotient Vn- The horizontal term yhor , 
as y varies at fixed s, sweeps out the horizontal subspace at s and these subspaces, 
as s varies, forms the horizontal distribution associated to a connection on the 
principal 5 ll -bundle S cm —► CP^ 2 . It is a basic fact about principal G-bundles 
with connection that the union of the horizontal spaces for the connection forms 
a G-equivariant vector bundle over the total space, and the quotient of this vector 
bundle by G is canonically isomorphic to the tangent space to the base space. 
Writing [s,j/] to denote the S 1 -equivalence class of the pair (s,y) we see that the 
set of all [s,2/hor]’s forms TCP W ~ 2 . Now s, together with (y? lor ,u, J) determine 
y uniquely. It follows that the map [s,j/] ([s, yhor], (v, J)) is a vector bundle 

isomorphism between the vector bundles ( S cm x E C m)/S 1 and TCP^ -2 x R 2 over 
CP^ -2 . The radial scaling coordinate r “goes along for the ride” without change. 
QED 


Because the decompositions of equations (23 24) are orthogonal and the second 
decomposition is scale invariant it follows that total kinetic energy decomposes as 


(25) 


r-/ \ 1 U IT Kshape([s : yhor]) 

A M)= “+2UT +-7- 


i y j 1 

i't + t 


K. 


shape 


The final term K s h a pe is formed by computing the squared length of the horizontal 
factor yhor and is canonically identified with the kinetic energy of the standard 
(Fubini-Study) metric on the shape space CP^ -2 . 

Remark. The kinetic energy decomposition (25) shows that for J / 0 the 
manifolds M mt (H 0 , J) is already closed in Vn so that 


( 26 ) 


M{h,J) = M int {h,J) 
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Indeed, the energy equation rh = H shows that U > \J 2 /r + 0(r) holds on 
M lnt (h,J ) which shows that if for a sequence Pi £ M mt (H 0 , J) we have that 
r(pi ) —>• 0 then U(si ) -> oo so that the shape Si of these points pi are converging 
to the collision locus C C CP^ 2 on the shape space. But we deleted C in forming 
Vn- 


5.1. Euler-Lagrange family in reduced coordinates. We follow Moeckel and 
look into what a planar central configuration family of section |4.3.1| such as the 
Euler or Lagrange family looks like in the coordinates of Vn- 

Let s cc be a planar central configuration and [s cc ] £ CP^ -2 the corresponding 
point in shape space. During the evolution of the associated family, this shape does 
not change. Only the size r and angle 9 of of the configuration changes. This size 
and angle change is specificied by A = A (t) = re 16 where A(f) solves the Kepler 
problem as per exercise [6j Since the shape does not change, the shape velocity yhor 
is identically zero along each of these solutions and so K s h = 0. Thus along such a 
solution 


r> 1 9 1 ?2 1 2 1 J 

A = -v 2 + - J 2 = -V 1 ^ - 

2 2 2 2 r 


(see eq. (25)) and the only variables which change are (r, u, J) among the full set 
of variables (r, [s, yhor], (u, J)) of Vn = [0, oo) x T(CP N ~ 2 \C)xR 2 But J = r 1 / 2 J 
and J is constant along solutions so the change of r and choice of J determines the 
change of J. So we can think of the only variables being is, r. 

Fix the energy h. We can then view the central configuration family as a one- 
parameter family of curves in the (is, r) plane, the parameter being the angular 
momentum J. Indeed the energy equation reads: 


rh = 


+ ^J 2 /r -U(s cc ). 


and since U(s cc ) is constant, this defines a one-parameter family of curves. We plot 
these curves in the is, r plane for various values of the angular momentum J below 
in figure |5T| 

Observe the rest point cycle in this picture: the closed curve passing thru the 
two equilibria. This curve is the union of two solution curves, a top arch which is an 
honest solution, and a bottom return curve. The top arch is the ejection-collision 
orbit first described when we described Lagrange’s solution: it explode out of total 
collision along the shape s cc achieves a maximum size and shrink back to triple 
collision. It connects the rest point s £ C having shape s cc and v = -\/2U ( s cc ) > 0 
the rest point s* having shape s cc and v = — yj2U(s cc ) < 0. This top arch lies on 
M(h , 0). The bottom ‘return road’ lies on the full collision manifold {r = 0, H = 0} 
and yields a return route from s* to s. This rest point cycle is the limit of the family 
of the periodic central configuration solutions with J ^ 0 as J — ?• 0. 

Notational Convenience. We have just used the symbol M(h, 0) C Vn for 
what used to be a submanifold of the phase space before quotient. We will continue 
to use the same notation for any G invariant submanifold or function on phase 
space before or after the quotient procedure. Thus we have: 


C,M 0 ,M(h),M int (h),M(h,J 0 ), etc. C V N - 
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Figure 6 . A central configuration family in u, r coordinates. The 
arch and ‘floor’ r = 0 comprise the rest cycle 


6. A GRADIENT (LIKE) FLOW! 

The dominant aspect of the flow on the full collision manifold Mo is that 
acts like a Liapanov function. 


Exercise 9. Use equations (13) to derive the identity 

1 


(27) 


v' = K- ^v 2 + H 


(See for example Moeckel |T5] , eq. (1.6).) 
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Exercise 10. Use the ‘‘Saari decomposition” of kinetic energy (eq j[25\ ) ) to show 
that 

K-^u 2 = K sh + ^J 2 . 

Conclude, using the previous exercise, that 

v = K sh + ^ J 2 > 0 on M 0 = {r = 0, H = 0}. 

You have proved much of 

Proposition 3. v' > 0 everywhere on the full collision manifold Mq. Moreover v 
is constant along a solution lying in Mq if and only if that solution is one of the 
equilibria. 

Remark. A flow is called “gradient-like” if it admits a continuous function / 
which is strictly monotone decreasing along all solution curves except equilibria. 
(See Robinson [24] p. 357.) The proposition thus asserts that the blown-up flow is 
gradient-like on the full collision manifold M 0 relative to the function / = —v. 

PROOF of PROPOSITION [3j In the exercise you proved that v' is positive every¬ 
where except at the points where it is zero. We must then show that any solution 
which lies on the locus v' = 0 is an equilibrium. We see that v' = 0 if and only 
if K sh = 0 = J. Now dJ/dr = —\vJ. (This holds both on and off the collision 
manifold). It follows that any solution starting on J = 0 remains on the locus 
J = 0. K s h([s, z/]) = 0 if and only if y^or — 0 in which case, both the yhor and the 
is term (from J = ( is,y )) in the decomposition of y are zero so that y = Xs with 
A € 1. Take inner products with s to get A = v. Now assume we have a solution 
curve (s(t),i/(t) lying on the locus J = 0,K s h = 0. Differentiating the equation 
y(r) = u(t)s(t) using the blow-up equations we see that y' = is's + vs'. But z/ = 0 
by assumption and s' = y — vs = 0 by the blow-up equations, so y' = 0 along the 
solution: our curve is an equilibrium. 

QED 

Pause for a moment. Reflect how different flow on the full collision locus is from 
a Hamiltonian flow on an energy level set. 

6.1. Moeckel’s manifold with corner into a manifold with a T. In |T5], at 
the beginning of section 2, Moeckel constructs a certain manifold with corners in 
preparation for perturbing the heteroclinic tangles lying on M(h , 0) into the realms 
of M(h,e). (He denotes his manifold with a corner by M 0 + and later simply M.) 
Dynamics on this manifold-with-corners is essential to our proof of theorem 1. I 
had a hard time making sense of this manifold. I rederive what Moeckel did in a 
slightly different way. I get a “manifold with a T” instead of Moeckels manifold 
with a corner. A “T ” is made out of two corners, or “L”s (one reflected relative 
to the other) joined along their vertical edge. One of these corners is Moeckel’s 
manifold with a corner and the other is a reflection of it. The vertex, or corner 
itself, is our good friend C the standard collision manifold. (Figure |6T] .) 

Recall that M(h) is a hypersurface in Vn, and as such is a manifold with bound¬ 
ary, whose boundary is our friend full collision manifold Mq = {r = 0, H = 0}. 

Definition 10. M(h) = M(h , 0) U M 0 C M(h). 
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M(h ) is a codimension 1 subvariety of the smooth manifold with boundary M(h). 
It is the zero locus of the function rJ restricted to M{h) and as such has two 
algebraic components : r = 0 which is our full collision manifold Mo, and J = 0 
which forms M(h, 0). The singular locus of M(h) is the intersection C = {r = 
0, J = 0} of these two components. All the rest point cycles described above 
associated to the central configurations lie on this M(h). M[h) is comprised of two 
“manifolds with corners”, namely {rJ = 0, J > 0} and {rJ = 0, J < 0}. The first 
of these is Moeckel’s. 

M(h ) is to be viewed as the limit as J —>■ 0 of the manifolds M(h , J). 



c 


<--0 


Figure 7. M(h) 


inside M{h ) is the zero level set of rJ. 


Proposition 4. For Sc8 a subset of the line of angular momentum values, set 
M int (h, S ) = U JeS M mt (/i, J). Then M(h) = n e>0 M int (h , (-e, e)). 
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The proof of the proposition follows in a routine way from our expressions for 
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energy, from rh = H, J = r~ 1 ^ 2 J and the kinetic energy decomposition of eq 
It is useful to recall, eq (261 that the M(h , J) = M mt {h , J) are closed for J ^ 0. 

As an alternative to the description of the proposition, we can either let J -A 0 
from above or below. Set 


and 


M + (h) = lim M(h,J) 

J - >o+ 


M_Vh)= lim M(h,J). 

J —y 0 


Then one can show without difficulty that 

M(h) = M + (h) U 

with M+(h) ={j)£ M(h), J > 0} and M-(h) = {p £ M(h), J < 0} being the two 
manifolds with corners described earlier, Moeckel’s manifold with a corner being 
M+. 

What is a manifold with a ‘T’? Suppose we have two real-valued functions x, y on 
an n-dimensional manifold Q such that 0 is a regular value for both functions and 
(0, 0) is a regular value of the map (; x , y) : Q — > R 2 . Then the locus {xy = 0, y > 0} 
is a manifold with a T. Its singular locus is {x = y = 0}. A manifold with a T is 
locally diffeomorphic to the prod uct of the “upside down T” xy = 0,y > 0 in the 
xy plane, by an R n_2 . See figure 6.1 


6.2. Finishing up the proof of theorem 1. The idea of Moeckel is that hy¬ 
perbolic structures persist on perturbation, and that the various stable-unstable 
connections between Euler and Lagrange central configuration points on on M(h ) 
are sufficiently “hyperbolic” that they persist into M(h , e) for e yf 0 small. Nonzero 
angular momentum is needed to get orbits connecting from R's to R* in finite time 
since the rest cycle of figure [1 lT| takes infinite blown-up time. Moeckel cannot carry 
out the ‘perturbation of hyerbolic” idea literally because he cannot establish the 
needed hyperbolicity or transversality. Instead, following an earlier idea of Easton, 
he replaces hyperbolicity by a weaker notion of “topologically transverse” between 
collections of “windows” transverse to the flow. This notion is sufficiently flexible 
and stable to allow Moeckel to perturb the various formal connections to get actual 
orbits realizing walks in the abstract graph introduced in section [3] By following 
the details of his proof, three decades later, we were able to verify that his realiz¬ 
ing solutions when projected onto the shape space do indeed stay C 0 -close to the 
concrete connection graph as described in section [3] 

The hypothesis of equal or near equal masses is needed to insure that (some 
of) the eigenvalues for the linearization at the Euler equilibria are complex. This 
complexity implies a “spiralling” of the Lagrange stable/unstable manifolds around 
the Euler unstable/stable manifolds and is needed to insure that all connections in 
the abstract connection graph are realized. 
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